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Abstract. We investigate the viscoelastic properties of entangled networks of semiflexible polymers. At 
intermediate time scales the elastic response of these networks to shear deformation is described by the 
plateau modulus G. Different scaling laws with polymer concentration c have been proposed based on the 
assumption that the deformation field is affine on all length scales. We develop a numerical approach that 
allows to calculate the modulus via free energy changes for both affine and non-affine deformations. The 
non-affine deformation field is obtained by a free energy minimization. Our findings allow for a confirmation 

of a power law G oc c 7 ^ 5 ^ 1 ^ 5 with polymer concentration c and persistence length l p and furthermore 
quantify the systematic deviations due to the affinity assumption. 



PACS. 87.16.Ka Filaments, microtubules, their networks, and supramolecular assemblies - 
Properties - 62.25.-g Mechanical properties of nanoscale systems 



87.19.rd Elastic 



1 Introduction 

Semiflexible polymers have always been an interesting test- 
ing ground for concepts of statistical physics, since ex- 
citable energies are on the order of k^T, fluctuations are 
important and entropic effects compete with bending en- 
ergy. However, with a few exception like the tangent-tangent 
correlation function [IJ[2] , moments of the end-to-end dis- 
tance [3] in the wormlike-chain model [US] or the end- 
to-end distribution function [5] in the weakly-bending rod 
approximation, a multitude of interesting single polymer 
properties are not accesible to analytic solutions. This ap- 
plies even more so to many-polymer systems like networks 
or solutions. In these cases, inter-polymer interactions or 
topological constraints complicate the theoretical treat- 
ment and the mere number of constituents and degrees of 
freedom clearly renders any analytical attempt to calcu- 
late the partion sum unfeasible. Our line of approach is 
therefore to reduce the many-body problem to a single- 
polymer description. A prominent example of such an ap- 
proach is the tube model introduced by de Gennes !6J and 
Doi and Edwards PQ that describes the combined effect 
of neighboring polymers by a mean-field potential. These 
theoretical models for single polymers at hand, one can 
proceed to describe macroscopic properties of a polymeric 
material. Here it is essential to build on the microscopic 
constituents without loosing emerging collective proper- 
ties of the macroscopic material. 
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These material properties of polymer networks are of a 
stunning variety and complexity. Spurred by the biologi- 
cal importance of networks of filamentous actin (F-actin) , 
intensive experimental research has recently focussed on 
this model system. In the presence of cross-links polymer 
solution form permanent gels or bundles with different 
elastic properties determined by a competition of bending 
and stretching modes jTlfHlfffifTO] . Recently, also the im- 
portance of non-affine deformations has been pointed out 
[llj . For purely entangled solutions a strong dependence 
of both the storage and the loss modulus on frequency 
was observed [12,13,14.15]. Furthermore, similarities to 
glassy systems were reported [T6lfT7] . Upon application 
of larger stress a non-linear regime was investigated and 
shear stiffening of the network was observed pl8l[T9] . The- 
oretical analysis identified different scaling regimes for the 
moduli with frequency [5D] . Regarding the scaling of the 
plateau modulus with concentration Isambert and Maggs 
derived a power-law of 7/5 from a simple scaling argu- 
ment |21j based on the deformation of confinement tubes. 
Other theories attribute the material response to the sup- 
pression of undulations by stretching [Tll2"2"] . They predict 
a scaling with a considerably larger exponent, but seem 
to disagree with experimental data 23,24.25 . Recent ex- 
perimental work [26 claims to have verified an exponent 
of 4/3 predicted by an elastic medium theory [14]. Com- 
mon to all these approaches is the assumption that the 
macroscopic deformation field is assumed to be affinely 
transmitted to all length scales. We will present an ap- 
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proach that permits to go beyond this assumption and 
investigate resulting differences. 

We proceed as follows: in Section [2] we define the sys- 
tem under consideration. We introduce the Hamiltonian 
of the network, review the simplifications that lead to the 
tube model and prior work on the plateau modulus, and 
emphasize the mean- field nature of the tube model and the 
assumption of affine displacement. In Section [3] we present 
our approach to numerically compute the free energy of 
the system by a reduction to two dimensions. We explain 
the effect of global shear on the microscopic constituents 
of the network and introduce a free energy minimization 
procedure that results in a non-affine deformation field. 
We proceed in Section [4] with the presentation and inter- 
pretation of our results before we conclude in Section [5] 



2 System Definition 

We consider a network of semiflexible polymers that only 
interact via a hard-core potential and thus constrain each 
other topologically. The polymer density is given by the 
number v of polymers of length L per unit volume. The 
stiffness k of the polymers gives rise to a persistence length 
of l p = n/k B T. The configuration of the i-th polymer in 
space, i"i(s), is parameterized by arc length s and the av- 
erage distance between the networks constituents is char- 
acterized by the mesh size, £ = y/3/vL. Concerning the 
mechanical response of the network, we are interested in 
the dynamic processes that occur after a deformation has 
driven the system out of equilibrium. These relaxation 
processes occur on different time scales. While generally 
every stress can relax by reptation of the polymers, this 
process is dramatically slowed down due to topological 
constraints in crowded environments 27,28 . On interme- 
diate time scales relevant for the plateau modulus it can 
be assumed that the constraints imposed by surround- 
ing polymers can not be overcome and that the center 
of mass of all filaments does not change substantially. A 
given polymer i is then described by the worm-like chain 
model [H[3] with a Hamiltonian that has contributions 
from intra-polymer bending and from interactions with 
the neighboring filaments in the solution. This can be writ- 
ten as 
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ds 



ds 2 
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(1) 



where the function 0fj formally describes the hard-core 
interaction between the polymers i and j under the initial 
topology I. Note that this function does not only depend 
on the configurations r,;(s) and ij-(s) of the polymers like 
a conventional hard-core potential that would always be 
zero in a system of mathematical lines without excluded 
volume. Instead it crucially depends on the inital topol- 
ogy and rules out polymer crossing by returning an infinite 
energy if the two polymers interpenetrate Q. The Hamil- 
tonian of the complete system is obtained as H = J^i Hi ■ 

For times far larger than the relevant time scale for the 



Given the form of the Hamiltonian |(TJ) a calculation 
of the free energy F from the partition sum Z as F = 
~k B T\nZ is obviously not feasible as the partition sum 



Z = II l V[r t (s)}exp[-H/(k B T)} 



(2) 



amounts to multiple path integrals of a highly convoluted 
integrand. 

A simplifying description of the system was proposed 
with the famous tube model 6,1 where the combined 
effect of the fluctuating neighbor polymers on a single test 
polymer is described by an effective harmonic potential. 
While the tube model is the foundation for theories for 
different properties of semiflexible polymer networks like 
tube diameter [2911141150] or viscoelasticity [21] , it has to be 
kept in mind that it only provides a mean-field description 
of the microscopic constituents. In the remainder of this 
section we will review how free energies and mechanical 
properties can be derived from the tube model and point 
out possible shortcomings of this coarse-grained frame of 
description. 

The tube model can be applied to solutions of semi- 
flexible polymers were the confinement of a single poly- 
mer by its neighbors is sufficiently strong to guarantee 
that the transversal undulations of the polymer do not 
deviate far from an average contour in space - the tube 
backbone r°(s). This is the case, if persistence length and 
polymer contour length are substantially larger than the 
typical void spaces in the mesh of surrounding polymers, 
thus L, l p ^> £ as e.g. given for most F-actin networks. 
The complex sum in the second term of the Hamiltonian 
(JlJ can then conveniently be substituted by a harmonic 
potential with average strength 7 and minimum at the 
tube backbone: 



ds 



d 2 r(s) 
ds 2 



+ l(r( S )-r (s)) 2 



(3) 

As pointed out by Odijk 31j it is instructive to intro- 
duce an additional length scale Ld ~ {vL)~ 2 t b l]J b known 
as deflection or Odijk length. While the length scales L 
and lp describe the properties of the single polymer and 
the length scale £ describes the network, the deflection 
length Ld captures the interaction between both. It can 
be interpreted as a measure for the distance between two 
collisions between the encaged polymer and the tube walls 
and therefore the number of collisions of a polymer is given 
as L/Ld- This is also reflected in the free energy cost AF 
that arises from the restriction of the test polymer to a 
tube and is obtained by a path integration of ((3]) over all 
polymer configurations ^2] as 



AF = V2k B T— 
Ld 



(4) 



contraints can be overcome by reptation. This signifies that 
long time averages will reproduce the results obtained from an 
ensemble average with respect to the free polymer Hamilto- 



plateau modulus the contributions Ofj vanish as topological nian. 
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This signifies that at a spacing of Ld between two collisions 
every of the L/Ld contact points between polymer and 
tube contributes one k^T to the confinement free energy. 

Having derived the free energy in the coarse-grained 
tube model, the next step is to analyze the change in free 
energy at mechanical deformation to obtain the plateau 
modulus. As the polymers are described in terms of their 
tubes, it is obvious to investigate the effect of deformation 
on the tubes for which the free energy is known as rea- 
soned above. Together with a scaling law d oc c~ 3 ^ 5 l p ^ 5 
for the tube diameter d derived by Semenov [29] this line 
of reasoning was first used by Isambert and Maggs [21] to 
establish a scaling relation between plateau modulus and 
concentration. They argue that the macroscopic shear de- 
formation is affinely passed down to the tubes that are 
compressed or stretched depending on their orientation 
to the shear. The resulting change of the tube diameter 
causes a change in the deflection length Ld and with the 
help of ((U) the resulting modulus scales as 

G oc c 7 / 5 ^ 175 • ( 5 ) 
The same scaling was also obtained by other descriptions, 
e.g. by a modified Onsager theory for confinement tubes 
|23] . In an endevour to arrive at a quantitative theory for 
the plateau modulus Morse [14] has proposed two concep- 
tually different approaches: a detailed microscopic descrip- 
tion of the topological constraint imposed by neighboring 
polymers leads him to the prediction of a modulus scaling 
with c 7 / 5 and a quite different approach yields a scaling 
of G oc c 4 / 3 derived by a self-consistent treatment of the 
network as an elastic continuum. Since these values are 
numerically quite close, a decisive distinction between the 
models has not yet been possible with the accuracy of 
available experimental data. 

It has to be kept in mind that all these theoretical 
approaches are implicitly build on the assumption that the 
tube contour deforms affinely with the macroscopic strain. 
This is evidently only a very coarse-grained description 
of the system's response. While stress relaxation by slow 
processes like reptation is obviously not relevant on the 
time scale of the plateau modulus, it is however possible 
that faster relaxation processes cause a tube contour that 
differs from the contour obtained by affine displacement. 
Our goal is to implement this relaxation processes by a free 
energy minimization in a numerical solution of the plateau 
modulus and investigate the quantitative and qualitative 
differences to the affine model. The detailed setup of this 
approach is discussed in the following section. 

3 Numerical Solution 

Our approach is to obtain a numerical solution of the par- 
tition sum ([2]) for a test polymer in a typical network of 
semiflexible polymers by averaging over all allowed config- 
urations of neighboring polymers. The advantage of this 
approach is a microscopic description of the Hamiltonian. 
In contrast to the tube-model that only provides a coarse- 
grained description of the surrounding polymers, it ac- 
counts for the detailed interactions in a given realization 



of disorder. This permits to investigate the effect of local 
non-affine deformations of the encaged test polymer. Since 
the distribution of obstacle polymers around a given test 
filament is quite heterogeneous, it is expected that these 
non-affine deformations result in a lower global free en- 
ergy. Our aim is to find this free energy minimum by a 
numeric minimization procedure. 

3.1 Reduction to 2D 

To reach this goal, we start by decomposing the transverse 
undulations of the test polymer into two independent com- 
ponents as previously described [33j . For one component 
the Hamiltonian thus simplifies to the description of a 
two-dimensional polymer in a plane surrounded by a cer- 
tain number iVobs fluctuating point- like obstacles (see Fig. 
[1]). The point-like obstacles are subjected to a harmonic 




Fig. 1. Fixed polymer in an array of fluctuating point obsta- 
cles (black points). The interaction between the polymer and 
the obstacles as the fluctuations of the j-th point obstacle are 
hindered by the polymer at a distance d" (see inset). 



potential with strength 7 around an equilibrium position 
p° with j = 1, .., The parameters N ohs and 7 can be 
chosen to self-consistently represent a network of a specific 
concentration |30] . Therefore the system is completely de- 
scribed by the two-dimensional contour r(s) of the test 
polymer and N b s two-dimensional vectors pj describing 
the positions of the obstacles in the plane. The Hamilto- 
nian thus reads H = H p + J2f=T H ° hs where RV is the 
bending energy contribution from the polymer 

and the H° hs are the contributions from the obstacle points 
H? s = \ (p, - P, ) 2 + 0[r( S ),p,, Pj o ] , (7) 
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where the topological constraint of uncrossability is again 
described by a function O that returns an infinite energy 
if the polymer and a point obstacle cross. In contrast to 
the three-dimensional case Eq. ([1]) this function can now 
easily be grasped geometrically as we will show below. In 
calculating the corresponding partition sum we have to 
solve 

„ „ N abB N ooa 

Z= V[v{s)] J J] d V] eM-pH P ]zM-PY. H ? S } ■ 

(8) 

While the description has now reduced to a single path 
integral, an analytical solution is still complicated by the 
topological constraints. However, for a specific polymer 
contour, the integration over the degrees of freedom of 
the obstacles is straightforward and can easily be carried 
out analytically as the only topological restriction for each 
point obstacle is posed by the test polymer. Since the poly- 
mer is mostly straight on the length scale of the typical 
fluctuation width of an obstacle, we assume that integra- 
tion over the obstacle potential is performed only in the 
half-space that is limited by the test polymer at a nearest 
normal distance rf™ as depicted in Fig. [TJ The partition 
sum is then written as: 



, N obs 

Z = I V[r{s)] exp[-/3H p ] ]J -erfc 
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(9) 



To solve the remaining path integration we chose to apply 
a saddle point approximation in which we first assume the 
test polymer to be immobile. Then we find the fixed con- 
tour that maximizes the partition sum and thereby mini- 
mizes the free energy and finally add fluctuations around 
this minimum. In the first step we are faced with the mini- 
mization problem depicted in Fig.[JJ An immobile polymer 
with associated bending stiffness is placed in an array of 
fluctuating obstacles. The free energy is composed of the 
bending energy of the polymer and the entropic contribu- 
tions from the obstacles. It is obtained from the partition 
sum from Eq. and is a function of the polymer contour 
r(s) alone: 



N obs , 

F(r(s)) = H p (r(s)) + ^ -k B T\n ( -erfc 
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(10) 

We are now looking for the contour r°(s) of the polymer 
that minimizes the free energy for a given initial setup of 
obstacles. 



3.2 Free Energy Minimization 

Technically, this contour is obtained as follows. We start 
from a given initial polymer configuration and choose N m i n 
nodes as a discretization along its contour. This reduces 
the required minimization to N m i n dimensions and per- 
mits to obtain a feasible computation time by a suitable 
choice of N m i n . To find the minimum of the free energy, 



Eq. (|10p , we move these nodes certain distances transverse 
to the present polymer contour, and define the new trial 
contour as a cubic spline through the new node positions; 
here one has to make sure that total contour length is kept 
constant. With this new contour both the bending energy 
and the entropic contribution to the free energy is com- 
puted according to Eq. (fT0|) . These steps of moves trans- 
verse to the immediately preceding contour are repeated 
until the minimum of the free energy is reached. The min- 
imization algorithm we used is based on the AMOEBA 
[33] implementation of the Nealder-Mead method 35J. 



3.3 Mode Representation 

Now that we have found the polymer contour of lowest free 
energy, we proceed to add the transversal fluctuations of 
the polymer in form bending modes. The contour of min- 




Fig. 2. Modes k — (red) and k = 3 (green) around the 
contour of minimal free energy (black). 



imal free energy r°(s) can be interpreted as the backbone 
of the test polymer's confinement tube or the contour with 
the highest probability. All deviations from this contour 
have higher free energy and thus a smaller probability. We 
model the thermal undulations of the polymer around this 
tube backbone by cosine modes ttfc(s) in the form 



r ±( s ) = ^2 u k(s) = ^ A fc COs(^-) 



(11) 



where k is the mode number and is the mode am- 
plitude. In this representation the mode k = is simply 
a transversal displacement of every point of the polymer 
normally to r° . A visualization of this mode and the mode 
Mfc=3 is depicted in Fig. [2] For a specific mode we can mon- 
itor the resulting free energy as a function of the mode 
amplitude A% as exemplary shown in Fig. [3J The result 
is a harmonic function F(Af.) = uj^/2A^. where u>u can 
be determined from the plots for every mode. The Hamil- 
tonian for a certain contour that is in the representation 
(fTTj) fully characterized by the set of coefficients {Af.} is 
then given as 



H({A k }) 



(12) 
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Fig. 3. Increase of free energy with amplitude for different 
modes in a given random array of obstacles. 



the deformation field follows the macroscopic stress on all 
length scales. Since at the time scale of the plateau mod- 
ulus the tube is the relevant quantity, the tube centers 
or backbones are displaced affinely with the global shear. 
Translated to our two-dimensional plane of observation, 
this signifies that the tube contour r°(s) and the centers 
of the obstacle tube p® are deformed affinely as depicted 
in Fig. [¥] (top). The free energy of this new configuration 
can be calculated as shown above. The modulus Ga to be 
determined from this resulting free energy change should 
be equivalent to the modulus obtained from any theoret- 
ical treatment that is based on the assumption of affine 
displacement. 



The desired partition function can thus easily be obtained 
by integration over the mode amplitudes 



Z = Y[dA keX p[-0u k Al/2]=ll^-. (13) 

The free energy cost of confinement is obtained by calcu- 
lating for every mode the free energy difference between a 
confined and a free polymer. Since lou oc fc 4 (compare Fig. 
[3]) this difference becomes small quickly with increasing 
mode number k such that it suffices to compute the free 
energy difference for the first few modes. 

3.4 Shear Deformation 

Having developed an approach to calculate the free en- 
ergy of a fluctuating polymer in an array of fluctuating 
topological constraints, we can proceed to investigate the 
free energy change as a reaction to shear deformations. 
A simple example of such a deformation is a global shear 
deformation of a macroscopic sample. If the sample is an 
equilibrated network of semiflexible polymers, the result 
of the shear deformation will be a rise in free energy and 
consequently a force counter-acting the deformation. The 
system is thus perturbed by the deformation and brought 
to a non-equilibrium state which will immediately be fol- 
lowed by relaxation processes. These relaxation processes 
occur on very different time scales for the different length 
scales in the network. This is the reason for the frequency 
dependence of the modulus. 

For very long time for instance, the network is able 
to completely relax the deformation stress by reptation 
thereby recovering the equilibrium value of free energy and 
resulting in a vanishing modulus. We can assume that at 
the time scale of the plateau modulus the encaged poly- 
mers have completely experienced their immediate sur- 
roundings but no large scale relaxation by network rear- 
rangement has occurred. Thus in measuring the plateau 
modulus the tube has sufficient time to form before the de- 
formation field changes again. This argument is the foun- 
dation for the assumption of affine displacement of the 
tube's contour and size. As it is unknown how exactly the 
macroscopic stress is passed on to the microscopic con- 
stituents of the network, it is commonly assumed that 




Fig. 4. Different levels of affinity in shear deformations, 
(top) Obstacle fluctuation centers and tube backbone follow 
the macroscopic shear deformation (large red arrow) affinely. 
(middle) While the obstacle fluctuation centers follow the 
macroscopic shear deformation affinely, the tube backbone can 
deviate from the affine deformation field in order to minimize 
the global free energy, (bottom) Also the obstacle point can 
react by non-affine deformation to the macroscopic shear. 

If we take a closer look at the processes at tube for- 
mation, it becomes clear that the tube obtained by affine 
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deformation is not necessarily a valid description of the 
physical reality. To this end we relax the assumption of 
afhnc displacement in such a way, that we now only let the 
obstacle fluctuation center p® deform affmely as illustrated 
in Fig.[4](cenier). If this new configuration of the obstacle 
array is taken for granted, we have to ask the question how 
the test polymer encaged by the obstacles reacts to this 
conformational change. Out of all possible deformations of 
the test polymer only the one with the lowest free energy 
will actually be realized. Obviously, this new tube con- 
tour will only in very few cases equal the tube contour ob- 
tained by affine displacement of the original contour. The 
resulting free energy of this deformation will thus be less 
or equal to the free energy obtained by affine deformation 
and we can therefore also state that the resulting modulus 
Gna < Ga- Technically, the free energy difference of this 
non-affine deformation of the tube contour is obtained by 
displacing the p°- affmely with the macroscopic stress (see 
[A| and then applying again the free energy minimization 
as explained above to find the new tube contour. 

Of course, also the fact that we deform the obstacle 
points affmely implies an assumption. In the actual phys- 
ical system the obstacle points and with it the tube cen- 
ters of the neighboring polymers are free to change their 
position in order to reach a global state of lower free en- 
ergy (see Fig. [4] [bottom). As these neighboring polymers 
however couple to other polymers outside our plane of ob- 
servation, the incorporation of this feature would be tan- 
tamount to a minimization in all degrees of freedom of 
the network. This is obviously out of range of a numerical 
solution. The resulting modulus of a complete free energy 
minimization is the modulus G that would be observed in 
experiments. The modulus determined by our approach 
constitutes an upper bound for the experimental values 
and thus G < Gna- 



4 Results 

In the previous section we presented an approximative nu- 
merical solution to the problem of finding the free energy 
change under shear deformation of a single probe poly- 
mer. This is obviously a quantity that can not be observed 
experimentally, but it can serve as the basis for the calcu- 
lation of the macroscopic plateau modulus. To obtain this 
observable we have to add up the contributions from all 
polymers in the network under consideration. One single 
specific polymer is described in terms of the two dimen- 
sional plane of observation whose orientation in space is 
described by a set of angles (6, <fi, tp) as described in Ap- 
pendix [A] Consequently, we have to perform an average 
over these isotropically distributed angles and furthermore 
we have to average over the quenched disorder that is gen- 
erated by the different configurations of point-like obsta- 
cles in the observation plane. With the shear parameter 
r this procedure finally returns an average free energy 
function AF(r) — gT 2 /2 from which the macroscopic 
modulus is obtained as G — 2vg. The factor 2 stems from 
the fact that every polymer is described by two planes 



of observation corresponding to the two components of 
transverse fluctuation. 

For a single polymer in a plane and one specific realiza- 
tion of obstacle disorder the resulting free energy function 
is exemplary shown in Fig. [5] (top). Since the absolute 
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Fig. 5. Free energy change with shear r for three different 
realizations of a test polymer in a network. 



value of the free energy differs strongly with the actual 
obstacle configuration all plots have been rescaled to the 
free energy value at r = 0. Obviously, for a single poly- 
mer in different specific realizations of obstacle disorder 
and different orientations of the plane of observation to 
the applied shear the resulting form of the free energy 
is highly variable. Furthermore, the free energy minimum 
is in general not at the point of zero shear. This feature 
however, should of course be fulfilled for the free energy 
function that is obtained by summing up all constituents 
in an macroscopic sample at equilibrium. We chose to use 
this requirement as a verification for an sufficient sampling 
over disorder. The location of the accumulated free energy 
minimum initially strongly oscillates with the number of 
samples but finally converges to P = 0. For every data 
point we average over a sufficient number of disorder sam- 
ples until this criterion is fulfilled. 

An example of the resulting averaged free energy func- 
tion is shown in Fig. [5] As expected the minimum is at 
zero shear where the system is at equilibrium. At the ap- 
plication of small shear the free energy rises in a harmonic 
fashion which would entail a linear restoring force in an 
experimental measurement. Furthermore, it can be seen 
that the free energy obtained in the affine approximation 
is always above the non-affine free energy that was ob- 
tained by the minimization procedure explained above. 
We determine the modulus by an harmonic fit at small 
shear strains. At higher strains however, the free energy 
function is no longer faithfully described by this fit, but 
features a stronger slope. This signifies the onset of non- 
linear forces. 
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Fig. 6. Free energy change obtained by averaging over 
quenched disorder and orientation where the lower free energy 
curve is due to energy minimized non-affine deformations. Har- 
monic fits are only a valid approximation in the linear regime 
at small shears. 

4.1 Affine vs. Non-Affine 

We determined the resulting plateau moduli as a function 
of different system parameters. Fig. [7] illustrates the scal- 
ing of the modulus with polymer concentration c. Both the 



tion that is obtained by permitting the encaged polymer 
to find its tube of minimal free energy leads to a lower 
modulus. Comparing the moduli from our affine calcula- 
tion with the prediction for the absolute plateau modulus 
by Morse's "Binary Collision Approximation" [TS] . shows 
sound consistency. In the realm of the restriction of affine 
displacement, our work can be seen as a numerical con- 
firmation. However, the moduli obtained by experimen- 
tal measurements [2 3) are considerably lower and prove 
that the physical reality is closer to a non-affine deforma- 
tion field. Of course, the detailed nature of this field is 
not accessible to experiments but our data suggests that 
the proposed model of an affine displacement of neighbors 
combined with an non-affine displacement of the tube is 
an appropriate approximation. The non-affine moduli ob- 
tained by this approach only show slight overestimations 
of the experimental results and it can be argued that this 
is due to possible additional non-affinities in the obstacle 
displacement. 



4.2 Scaling with Persistence 

Finally, we determined the scaling of the plateau modu- 
lus with persistence length l p . The decrease of the mod- 
ulus with increasing polymer stiffness is depicted in Fig. 
[8] and shows good agreement with a power law of —1/5. 
The persistence dependence represents a sensible method 
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Fig. 7. Moduli resulting from affine and non-affine displace- 
ment of the tube contour as a function of actin concentra- 
tion comply with a 7/5 power law. The affine modulus (filled 
squares) is in the range of the prediction by Morse [14] while 
the non-affine modulus (open squares) is a factor two to three 
lower but slightly above the experimental measurements (open 
triangles) by Hinner [23] . 



affine and the non-affine modulus show good agreement to 
a 7/5 power law with concentration. The non-affine modu- 
lus is considerably below the affine modulus. This confirms 
the initial assumption that a deformation field that as- 
sumes affine displacement on all length scales indeed over 
estimates the system's response. The non-affine deforma- 
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Fig. 8. The non-affine plateau modulus clearly shows a —1/5 
power law dependence on persistence length. 



to discriminate between competing models of viscoelastic- 
ity. Since the value of the concentration scaling exponent 
4/3 predicted by an effective medium approach [14] is nu- 
merically quite close to the exponent 7/5 predicted by 
most other theories, experimental accuracy does not al- 
low for a verification. The difference of the two concepts 
in the persistence length scaling exponents is considerably 
larger: —1/3 versus —1/5. Contrary to recent experiments 
[26j our data is clearly incompatible with an exponent oft 
— 1/5. We therefore conclude that the plateau modulus of 
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entangled networks of semiflexible polymers is correctly 
described by 

G oc c 7 / 5 /; 1 / 5 . (14) 

5 Conclusion 

We have presented numerical results for the plateau mod- 
ulus of entangled network of semiflexible polymers for a 
wide parameter range for both an affine and non-affinc 
shear deformation field. To this end we have developed 
an approach that permits to approximately calculate the 
partition sum and the free energy of a polymer network. 
This was achieved by analyzing the free energy of one 
component of the transversal fluctuations of a test poly- 
mer in a two dimensional reference frame. Averaging over 
disorder and all possible reference frames results in a mea- 
sure for the systems free energy. The approach allows to 
probe the system's free energy change and thereby me- 
chanical response to macroscopic and microscopic defor- 
mation fields. While existing theories for the modulus of 
polymer networks are based on the assumption of an affine 
shear deformation on all length scales, we challenged this 
assumptions. Indeed, it was found that the free energy 
of the affine deformation can be reduced by allowing the 
tube contour of a test polymer to minimize the global en- 
ergy. We observed non-affine moduli that agree well with 
experimental data, while existing theoretic predictions co- 
incide with the results of our considerably higher affine 
moduli. Furthermore, we clearly confirm a scaling of the 
plateau modulus with persistence length and concentra- 
tion as G oc c 7 / 5 lp 1 ^ 5 . Our results prove that shear de- 
formation of networks of entangled polymers has to be 
described in a non-affine picture and that affine theories 
systematically overestimate the mechanical response. The 
presented approach provides a numerical solution to eval- 
uate complex partion sums and has a wide applicability 
to rheology of polymer networks. Future applications can 
e.g. investigation of non-linear shear. 

We acknowledge support from the DFG through grant Fr 850/6- 
1, from the German Excellence Initiative via the NIM program 
and from the Elite Network of Bavaria through the NBT pro- 
gram. 



A Shear Deformation 

We work with two coordinate systems: a three dimensional 
real space system and a two dimensional system in the 
plane of observation, where the origin of both systems is 
one end of the initial test polymer. The orientation of the 
end-to-end vector R of an arbitrary test polymer in the 
three dimensional space is isotropically distributed. As de- 
picted in Fig. [5] it is described by the two angles 9 and 4>. 
To define a plane of fluctuations we need one additional 
angle ip. This plane is spanned by the vectors R and S 
with SIR. For tp = the vector S is obtained by ap- 
plying the same transformation to an vector parallel to 



the z-axis, that is needed to transform an vector paral- 
lel to the x-axis to R. Other values of ip are obtained by 
an rotation around the axis R. Graphically it is helpful 
to picture the observation plane as the plane that is ob- 
tained by applying two transformations to the x-z-plane: 
first a rotation around the z-axis by 4> and then a rotation 
around the axis R by ip. If R and S are normalized they 
correspond to the x and y-axis in the plane of observation. 

If we now apply a macroscopic shear deformation T(r) 
with shear parameter _T to the three dimensional system 
the obstacle points and the test polymer deform according 
to the action of the transformation T on their real space 
coordinates. In general this signifies that they leave the 
plane of observation spanned by R and S. It is however, 
self-evident that also the test polymer's fluctuations are 
subjected to the shear and therefore also the plane of ob- 
servation transforms according to T . This is tantamount 
to a transformation of the vectors R and S and leaves the 
transformed obstacles points in the new plane of observa- 
tion. We obtain the new plane coordinates in terms of the 
transformed unit vectors R' and S' that correspond again 
to x and y-axis. The former is obtained as R' = 7~(T)R 
and the latter is constructed as the component of T(r)S 
that is orthogonal to R'. 




y 

Fig. 9. Schematic illustration of the plane of observation. 
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